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■ The AA bond energy ASaa in AA hypernuclei is obtained from a G- 



matrix calculation which includes the coupling between the AA, 'EN and 
channels, as well as the effect of Pauli blocking to all orders. The Nijmegen 
NSC97e model is used as bare baryon-baryon interaction in the strangeness 



^ . S = —2 sector. The AA-HA^ coupling increases substantially the bond energy 

with respect to the uncoupled AA case. However, the additional incorporation 
of the channel, which couples simultaneously to AA and EN states, has 
a surprisingly drastic effect and reduces the bond energy down to a value 
closer to that obtained in an uncoupled calculation. We find that a complete 
treatment of Pauli blocking reduces the repulsive effect on the bond energy 
to about half of what was claimed before. 
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I. INTRODUCTION 

Double-strange A hypernuclei are nowadays the best systems to investigate the properties 
of the S — —2 baryon-baryon interaction. Emulsion experiments and subsequent analysis 
[1-4] have reported the formation of a few AA hypernuclei, aaHc, aaBc and ^B. Prom the 
resulting AA binding energies, a quite large AA bond energy of around 4-5 MeV emerged, 
contrary to expectations from SU(3) [5]. A series of theoretical works, based either on 
phenomenological S = —2 baryon-baryon interactions or realistic ones, have studied the 
properties of double-A hypernuclei for more than 30 years [6-23] . 

The recent finding at KEK of a new aaHc candidate having a AA bond energy of around 1 
MeV [24] has injected a renewed interest on this field. Unless new experiments for the other 
AA hypernuclei also give lower binding energies in the future, it is now an open question 
to reconcile theoretically the weak attraction found in %jjie with the stronger attraction in 
the other two heavier systems. Although some progress has been made in Ref. [25], where 
both short- and long-range correlations were simultaneously treated, further investigations 
are needed to completely settle this question. Fihkhin and Gal [26-28] report Faddeev- 
Yakubovsky calculations, complementary to those carried out in the earlier work of Ref. [16] 
but using the new Nijmegen interactions, not finding a simultaneous description of the ]^\Be 
and the new A^He binding energies. However, the nucleon Pauli blocking effect affecting, 
through the couphng to EN states, the AA interaction when the particles are embedded 
in the nuclear medium has not been considered in most of the earlier nor in these recent 
works [25-27]. As discussed in detail in Ref. [17], Pauli blocking reduces substantially 
the additional attraction to the AA binding energy induced by the AA — > EN conversion. 
Recently, an attempt to incorporate the Pauli suppression effect has been made in Ref. [29], 
where a second order Pauli correcting term is introduced in the intermediate states following 
the AA — > EN transition. The interaction used in that work is a two-channel (AA, EN) 
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Gaussian model, which imphcitly includes the EE coupling not only in the effective AA 
interaction but also in the AA E,N transition. In fact, the important role of the coupling 
to EE states has been recently pointed out in Ref. [30] and explicitly worked out for the 
Nijmegen interactions in the variational calculation of Ref. [31]. 

The purpose of the present work is to present a careful analysis of the role of coupled 
channels on the S = —2 baryon-baryon interaction in the medium, treating Pauli blocking 
effects to all orders in all possible transition channels. For practical purposes, most of the 
recent works [26-29] have used simple parameterizations of the new Nijmegen potentials 
in terms of the sum of a few gaussians. In contrast, we start from the original Nijmegen 
model NSC97e [5], as done also in Ref. [31]. In our approach, we solve the coupled-channel 
equation for the G-matrix in infinite nuclear matter, and derive from it the AA bond energy 
in finite nuclei. With respect to existing calculations our treatment of the finite system 
is very simple. This has the practical advantage of permitting us to explore in depth the 
different effects determining the AA bond energy, such as coupled channels or Pauh blocking 
to all orders. 



II. FORMALISM 

The AA bond energy AB\\ in AA hypernuclei is determined experimentally from the 
measurement of the binding energies of double- and single-A hypernuclei as 

ABaaH^Z) = 5aa(^aZ) - 2BA{i-'Z) . (1) 

A reasonable estimation of this quantity when rearrangement effects are small can be ob- 
tained from the value of the AA G— matrix element in a finite hypernucleus 

AB^aHj^Z) ^ -((0si/2)a(0si/2)a, J = I G I (0si/2)a(0si/2)a, J = 0) , (2) 

where the two A particles are assumed to be in the lowest single particle state of an appro- 
priate A-nucleus mean field potential. We will compute the above matrix element from the 
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infinite nuclear matter one in the following way. First, we construct the AA G— matrix in 
infinite matter by solving the well known Bethe-Goldstone equation which, in partial wave 
decomposition and using the quantum numbers of the relative and center-of-mass motion, 
reads 

{K q' L' S' J {A A)\G\KqLS J (AA)) = (Kq'L'S' J{AA)\V \KqLS J (AA)) 
+ E E / dq"q"^Kq'L'S'J{AA)\V\Kq"L"S"J{BiB2)) 

B1B2 L"S" 

^ Qb,bM,Q") 

X {Kq"L"S" J {BiB2)\G\KqLS J (AA)) , 

where the labels B1B2 run over AA, EN and EE intermediate states. The starting energy Q is 
taken equal to 2Ma — 2Ba{\^^Z) — ABaa{aa^) = 2Ma — i?AA(AA^), where the experimental 
value of Baa is taken for each hypernucleus, namely 7.25 MeV for aa^c [24], 17.7 MeV for 
\'^Be [1] and 27.5 for [3]. In this way, we are considering the interaction of each A particle 
not only with the nucleons in the nucleus but also with the other A particle. The nuclear 
matter density to be used in the Pauli operator Q is determined, for each hypernucleus, as 
the average nuclear density felt by the A particle in that hypernucleus. This is obtained by 
weighing the nuclear density at each point with the probability of finding the A particle: 

p = j p{r) I ^A{r) |2 d^r , (4) 

where p(r) is the nuclear density profile which is conveniently parameterized as 



being 



with a = 0.52 fm, R = I.UAn^^^ - OMAn fm and An the number of nucleons in the 
hypernucleus. The A wave function is obtained by solving the Schrodinger equation using 
a Woods-Saxon A-nucleus potential with parameters (Va, oa, Ra) adjusted to reproduce the 



experimental binding energy of the A in the single-A hypernucleus. For practical computa- 
tional purposes, from the resulting A r.m.s. radius we derive the oscillator parameter, 6a, 
of an equivalent harmonic oscillator wave function which will then be used in obtaining the 
finite hypernucleus two-body G-matrix elements of Eq. (2) from the nuclear matter ones 
displayed in Eq. (3). 

In the next step, we express the two-body ket state |(0si/2)a(0si/2)a, J = 0), built from 
the Osi/2 states of the equivalent harmonic oscillator potential, in terms of momentum and 
angular variables h,ji)A{k2, h^jijK-i J = 0) in the laboratory frame using 

|(0si/2)a(0si/2)a, J = 0) = j J dMA:2A;iA:^i?oo(&AA;i)i?oo(&AA:2)|(A;i,0,^)A(A;2,0,^)A, J = 0) , 

(7) 

where Rni{x) is the corresponding harmonic oscillator function. 

Finally, we express the two-body state with laboratory coordinates in terms of the states 
with variables in the relative and center-of-mass system, \KqLSJ{AA)), used in the solution 
of the Bethe-Goldstone equation 

|(^i,0,i)A(A;2,0,^)A, J = 0) = J dKK^J dqq^{KqOOO{AA)\kiQ^k2Q^, J = Q)\KqQQQ{AA)) , 

(8) 

where {KqOOO{AA)\ki0^k20^, J — 0) are the appropriate transformation coefficients [32,33] 
from the relative and center-of-mass frame to the laboratory system. We note that the only 
contribution comes from the partial wave ^5*0. Transforming the bra state ((0si/2)a(0si/2)a| 
in a similar way, one can finally evaluate the AA bond energy of Eq. (2) in terms of the 
infinite nuclear matter AA G— matrix elements. 

In Table I we summarize all the parameters that allow us to determine the relevant A and 
nuclear properties needed in the evaluation of the AA bond energy for the three hypernuclei 
studied in this work, %jjie, aa^® aaB- 
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III. RESULTS 



The diagonal ^Sq AA G- matrix element for zero center-of-mass momentum and zero 
relative momentum is shown in Fig. 1 as a function of the nuclear matter density for several 
starting energy values. As density increases, the G-matrix element loses attraction as a 
result of Pauh blocking which reduces the available phase space for the intermediate EN 
states. On the other hand, the G-matrix element gains attraction when the starting energy 
increases, since the coupling to intermediate states is then more efficient. We will return to 
this behavior when the results of finite hypernuclei are discussed. 

Table II displays our results for ABaa in AA^e, for various coupled-channel cases. The 
value of ASaa obtained from a calculation that neglects Pauli blocking effects, i.e., directly 
from the T-matrix, is also displayed with brackets. As expected, incorporating the coupling 
between the AA and EN, produces a drastic effect over the AA uncoupled situation, increas- 
ing AB\\ from 0.16 MeV to 0.78 MeV, a value that lies very close to the new experimental 
datum [24]. We note that, contrary to what it seems to be implied in Ref. [30], the couphng 
between the AA and EN channels is important even when the interaction is weak, as it is 
the case of the NSC97e potential used here which produces a scattering length of about 
-0.5 fm. Actually, in Refs. [17,30] the potential is adjusted for each coupled-channel case to 
reproduce a common value of the scattering length. Therefore, part of the coupling effect 
is embedded in the readjusted parameters. We now turn to analyzing the effect of the EE 
channel, located more than 150 MeV higher in energy from the AA and EN channels, and 
which has usually been neglected or taken in an effective way within single-channel (AA) or 
two-channel (AA, EN) interaction models. The results shown in Table II reveal, surprisingly, 
that the role of the EE channel is very important and reduces substantially the two-channel 
(AA, EN) value of A^aa down to 0.28 MeV, which is closer to the uncoupled single-channel 
result. Note that the repulsion found for the full coupled-channel G-matrix element around 
the AA threshold does not necessarily mean that the EE channel produces a more repulsive 
interaction. In fact, the (AA, EN, EE) Nijmegen model becomes so attractive that it even 
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supports a spurious deeply bound YY state around 1500 MeV below the AA threshold [31]. 
However, the size of the G-matrix will not be affected by the presence of this bound state 
since it lies very far away from the region of energies required by our model. The net effect 
around the AA threshold is that the full coupled-channel calculation has a smaller bond 
energy than the case in which only the AA and the EN channels are retained. 

Comparing the results of AB\\ with those between brackets, which have been obtained 
from a T-matrix calculation, one observes that Pauli blocking effects (non-existing in the 
single channel AA case) are quite important, especially when the three-channels (AA, EN 
and EE) are considered, reducing by half the value of AB\\. We note that our Pauli 
unblocked value of 0.54 MeV, obtained for the complete coupled channel calculation using the 
original Nijmegen potential NSC97e, is reasonably close to results obtained directly in finite 
hypernuclei but using effective Gaussian parameterizations fitted to the scattering length of 
the Nijmegen NSC97e interaction, namely ABaa = 0.58 MeV [26] and ABaa = 0.64 MeV 
[29] . We note that the recent variational calculation using the Nijmegen interactions quotes 
a slightly larger value of 0.81 MeV [31]. 

The results of Table II show that a proper treatment of Pauh blocking, neglected in 
most of the calculations using more sophisticated ways of treating the finite hypernucleus 
[16,25,26] , is needed to draw conclusions on the particular value of ABaa predicted by a given 
interaction. A first attempt to incorporate the Pauli suppression effect within the context of 
finite hypernuclei has been done recently in Ref . [29] , where a Pauli blocking term, correcting 
the phase space of intermediate EN states accessed via AA — > EN conversion up to second 
order in the effective interaction, is added. The A A bond energy is then reduced from 0.64 
MeV to 0.21 MeV, hence finding a Pauli suppression of 0.43 MeV, which is about twice the 
size of the reduction we find in the present work, namely (0.54 — 0.28) MeV = 0.26 MeV. 
The reason for the difference has to be found in higher order terms of the Pauli correction. 
Indeed, if, in the spirit of the procedure followed in Ref. [29], we truncate the series that 
defines the T-matrix in terms of the G-matrix, T = G + G{l/E — Q/E)T, up to second 
order in G, then the contribution of the Pauli blocking correcting term, G{l/E — Q/E)G, 
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amounts to 0.36 MeV. This is consistent with the value of 0.43 MeV quoted in Ref. [29] 
which was obtained with a shghtly modified effective interaction to fit the new ASaa value 
in aaHc. Moreover, we have checked that the series converges to our T-matrix result and, 
hence, to our complete Pauli correction of 0.26 MeV. The Pauh correction built directly in 
the finite nucleus in the full coupled-channel calculation of Ref. [31] is also small and of the 
order of 0.2 MeV. 

The scattering length for each coupled-channel situation is also shown in Table II to 
illustrate, as in other works [26,30], its correlation with the AA bond energy, which increases 
as the magnitude of the scattering length increases. This correlation is to be expected 
since oaa is proportional to the T-matrix and ASaa is proportional to the corresponding 
medium modified G-matrix. We would also like to point out that the scattering length 
changes substantially for each of the coupled-channel cases. This is apparently different 
from the results shown in Refs. [17,30] but, as mentioned before, in these later works the 
interaction is readjusted in each coupled-channel calculation to reproduce a common value 
of the scattering length. 

Finally, we collect in Table III the results of ASaa for the three observed AA hypernuclei. 
The role of coupled channels is qualitatively similar in the three hypernuclei: coupling the 
channel to the AA channel increases Ai?AA substantially, while the additional incorpo- 
ration of the EE channel reduces the binding also substantially, bringing the value of the AA 
bond energy closer to the uncoupled result. We also observe that the heavier the nucleus the 
smaller the binding, contrary to what one would be expecting from the present experimental 
results. The trend found here is a reflection of the behavior of the ^5*0 G-matrix element 
shown in Fig. 1. Inspecting the nuclear structure parameters for each hypernucleus shown 
in Table I, we see that the nuclear density for aa^g is the largest, slightly above 1.5po, hence 
this hypernucleus has the strongest Pauli repulsive effect. However, the starting energy D, is 
also the largest, which produces a gain in attraction. For the range of densities and starting 
energies explored by the three AA hypernuclear systems studied here, the dependence of the 
G-matrix on the starting energy is twice more important than that on the density. The net 
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effect is tfiat tfie largest AB\\ value is obtained for the lightest system. 

IV. CONCLUSIONS 

In this work we have obtained the bond energy AB\\ in several AA hypernuclei, following 
a microscopic approach based on a G- matrix calculation in nuclear matter using, as S — —2 
interaction, the recent parameterization NSC97e of the Nijmegen group. We have identified 
the AA bond energy with the ^5*0 AA G- matrix element calculated for values of the nuclear 
density and starting energy appropriate for each hypernucleus. 

Our simplified finite-nucleus treatment has allowed us to explore in depth the effect 
of the various coupled channels and the importance of Pauli blocking on the intermediate 
EN states, paying a special attention to the role of the EE channel usually neglected in 
the literature. Consistently with other works, we find that the coupling between the EN 
and AA channels has a drastic effect, increasing by about 0.6 MeV the calculated AB^ 
in %jjie with respect to a single-channel AA calculation. Surprisingly, the additional incor- 
poration of the EE channel yields a non-negligible reduction in the binding of 0.4 MeV. It 
would be interesting to explore the role of the coupling to EE states in other three-channel 
S — —2 interactions, such as the Nijmegen hard-core potential F [34]. Unfortunately, our 
momentum-space method can only handle soft-core interaction models. 

We have also explored, within the complete three-channel approach, the effect of Pauli 
blocking, which is often neglected or considered in a truncated way in previous works. With 
respect to a T-matrix calculation, our calculated value of ASaa in \/Jie gets reduced by 
0.26 MeV, about half of what was found on the basis of a second order Pauli corrected 
calculation [29]. 

Due to our simplified treatment of nuclear structure, we do not expect a quantitative 
agreement with experimental data for the three hypernuclei studied. However, from the bulk 
of studies of double- A hypernuclei available in the literature, it seems unreasonable to think 
that a proper finite nucleus calculation which incorporates consistently core-polarization 
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effects, migfit cfiange tfie calculated bond energies substantially enough to obtain a simul- 
taneous agreement with the data. In this respect, our results confirm, in accordance with 
recent cluster calculations [26], the incompatibility between the experimental binding ener- 
gies of the light double-A hypernuclear species. We note, however, that the disagreement 
would be reduced if the tt" weak decay of the AA^e ground state was assumed to occur to 
the first excited state of ^Be, as pointed out by Filikhin and Gal [26], hence reducing the 
bond energy in ]^/fie to about 1 MeV. A clarification of the experimental situation, through 
new experiments and analyses, is certainly needed in order to test the theoretical models 
and make progress in the field of doubly-strange systems. 
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TABLES 

TABLE I. Parameters of the Woods-Saxon A-nucleus 
potential (VA,aA, Ra), equivalent A oscillator parameter 
(6a), effective nuclear density (p) and G-matrix starting en- 
ergy (CI) for each AA hypernucleus. 





iAHe 


i°ABe 


13 r> 

AA^ 


Va [MeV] 


28 


28 


28 


OA [fm] 


0.59 


0.59 


0.59 


Ra [fm] 


1.60 


2.06 


2.67 


6a [fm] 


2.23 


1.89 


1.82 


p [fm-3] 


0.277 


0.181 


0.176 


n [MeV] 


2224.12 


2213.67 


2204.17 



TABLE II. AA scattering length and AA bond energy in 
^y^He, for various channel couplings. Results within brackets 
ignore Pauli blocking effects. 





UAA [fm] 


A^AA [McV] 


AA 


-0.25 


0.16 (0.16) 


AA,E:iV 


-0.84 


0.78 (1.02) 


AA, EN, SS 


-0.49 


0.28 (0.54) 
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TABLE III. AA bond energy in ^^He, j^^Be and X\B, 
for various channel couplings. Units are in MeV. 





iAHe 


10 T3 

AA^6 


13 R 

AA^ 


AA 


0.16 


0.0046 


0.11 


AA,EN 


0.78 


0.97 


0.96 


AA, EN, EE 


0.28 


0.22 


0.11 


EXP: 


i.oii°:i [24] 


4.2 ± 0.4 [1] 


4.8 ± 0.7 [3] 
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